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FINITE ELEMENT DIFFERENTIAL FORMS ON CUBICAL 

MESHES 

DOUGLAS N. ARNOLD AND GERARD AWANOU 

Abstract. We develop a family of finite element spaces of differential forms 
defined on cubical meshes in any number of dimensions. The family contains el- 
ements of all polynomial degrees and all form degrees. In two dimensions, these 
include the serendipity finite elements and the rectangular BDM elements. In 
three dimensions they include a recent generalization of the serendipity spaces, 
and new i7(curl) and H(div) finite element spaces. Spaces in the family can be 
combined to give finite element subcomplexes of the de Rham complex which 
satisfy the basic hypotheses of the finite element exterior calculus, and hence 
can be used for stable discretization of a variety of problems. The construction 
and properties of the spaces are established in a uniform manner using finite 
clement exterior calculus. 



1. Introduction 

In this paper we develop a family of finite clement spaces S r A k (Th) of differential 
forms, where 7/t is a mesh of cubes in n > 1 dimensions, r > 1 is the polynomial de- 
gree, and < k < n is the form degree. Thus, in 3 dimensions, the space S r A k (T n ) 
is a finite element subspace of the Hilbert space H 1 , -ff(curl), H(div), or L 2 , ac- 
cording to whether k = 0, 1, 2, or 3. For n — 1 or 2, the spaces were previously 
known, while in three (or more) dimensions, they are mostly new. Specifically, our 
construction yields a new family of .ff (curl) elements and a new family of H(div) 
elements on cubical meshes in three dimensions. Our treatment in an exterior 
calculus framework allows all the spaces and their properties to be developed to- 
gether. The spaces combine together in complexes satisfying the basic hypotheses 
of the finite element exterior calculus [5]. This means that, in addition to their 
use individually, they can be used in pairs, 5 r -|_iA (7h.) x 5 r A (7ft) in a variety 
of mixed finite element applications, with stability and convergence following from 
the abstract theory of [S] . Element diagrams for some of the spaces are shown in 
Figure [TJ The dime nsion of the shape function spaces are given in Theorem |3.6| 
below and are tabulated for n < 4 and r < 7 in Table [T] 

Since the development of the first stable mixed finite elements for the Poisson 
equation by Raviart and Thomas in 1977, such elements have proven to be power- 
ful tools for numerical computation. Their paper T1J introduces a family of finite 
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Legend: symbols represent the value or moment of the indicated quantities: 
• - scalar (1 DOF); ib - tangential vector field (2 DOFs); 
W - vector field (3 DOFs); - tangential component of vector field on edge 

or normal component on face (1 DOF). 



Figure 1. Element diagrams for S r A k in three dimensions for r < 3. 



element discretizations of the space H (div, fi) for a two-dimensional domain f2, one 
for each polynomial degree. Together with a corresponding discontinuous piecewise 
polynomial discretization of L 2 (fi), these RT spaces stably discretize the mixed 
variational formulation of the Poisson equation on fi. In [TT], versions of the RT el- 
ements were given both for meshes of by triangles and by rectangles. Besides the 
original application of mixed finite elements for the Poisson equation, the RT ele- 
ments can be used together with the standard Lagrange finite element discretization 
of iJ 1 (fi) to give a stable mixed finite element discretization of the vector Poisson 
equation curl curl u — graddivw = /, in which the vector variable u is sought in 
-ff(div) and approximated by RT elements, and the scalar variable a = — divu 
is sought in H 1 and approximated by Lagrange elements. The RT elements were 
generalized to 3-dimensions by Nedelec [9] , with separate generalizations giving dis- 
cretizations of H(div, f2) and of if (curl, f2), including both tetrahedral and cubic 
mesh generalizations of each. 
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Table 1. Dimension of S r A k (I n ). 



In [7], Brezzi, Douglas, and Marini introduced a second family of finite element 
discretizations of H(div, Q) in two dimensions, for both triangular and rectangu- 
lar meshes. These BDM element have also proven to be very useful. Nedelec [TU] 
generalized the BDM family to tetrahedral meshes in three dimensions, giving ana- 
logues both for iJ(div) and if(curl). He also defined H(div) and ii(curl) elements 
on cubic meshes in [10 . However these cannot really be considered analogues of the 
BDM elements, as they do not seem to lead to stable mixed finite element pairs. 
The generalization of the BDM elements to H(div) on three-dimensional domains 
(but not -ff(curl)), was also made by Brezzi, Douglas, Duran, and Fortin in 
The paper [B] also introduced an analogue of the rectangular BDM elements to 
cubic meshes in three dimensions. 

The finite element exterior calculus [4j [5] has greatly clarified the relation of 
many of these mixed finite element methods. In exterior calculus the space H(drv) 
is viewed as a space of differential forms of degree n — 1 in n dimensions, while 
if (curl) is a space of 1-forms, H 1 a space of 0-forms, and L 2 a space of n-forms. 
These spaces are connected via the de Rham complex, in which the fundamental 
differential operators grad, curl, div (and others in higher dimensions) are unified 
as the exterior derivative. We refer the reader to Table 2.2 in [4] for a summary of 



^The finite element discretization of ff(div) in [6] is identical with that of |10| in the case 
of tetrahedral meshes. However the interior degrees of freedom differ. This does not affect the 
stability analysis for the mixed Poisson equation, but only the ones given by 1 101 can be used to 
establish stability when discretizing a mixed formulation of the vector Laplacian. 
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correspondences between differential forms and vector fields. In [4j [5] two funda- 
mental families of finite element differential forms are denned on simplicial meshes 
in n dimensions. Tables 5.1 and 5.2 in [3] summarize the correspondences between 
these spaces of finite element differential forms and classical finite element spaces in 
two and three dimensions. The V~ A k family specializes to the Lagrange elements, 
the RT elements, and the fully discontinuous polynomial elements, for k = 0, 1, and 
2, respectively, in two dimensions, and to the Lagrange elements, Nedelec's gen- 
eralizations of the RT elements to H(cml) and to ff(div), an d the discontinuous 
elements for k = 0, . . . , 3 in three dimensions. Taken together, these spaces form a 
complex which is a finite element subcomplex of the de Rham complex: 

-> V~A° ^ 4 • • • 4 TV A" -> 0. 

The second family of finite elements discussed in [H [5] is the P r A k family. For 0- 
forms and n-forms this brings nothing new, just recapturing the Lagrange and fully 
discontinuous polynomial elements, but for < k < n, this is a different family of 
finite element spaces. For n = 2, k = 1 it gives the BDM triangular elements, and 
for n = 3, k = 1 and 2, it gives Nedelec's generalizations of these to H(div) and 
if (curl). They combine into a second finite element de Rham subcomplex 

-> P r A° -4 Vr-iA 1 -4 • • • 4 T r - n A n -> 0. 

Note that the degree decreases in this complex, in contrast to the preceeding one. 

We now turn to the construction of analogous spaces and complexes for cubical 
meshes. An analogue to the V~ A k complex of elements for cubical meshes may be 
easily constructed via a tensor product construction. For an explicit description in 
n dimensions and for all form degrees k, see [3]. This includes the tensor product 
Lagrange, or Q r , elements for 0- forms, the rectangular RT elements for 1-forms in 
2-D, and the 3-D generalizations of them given in [5]. The space for n-forms (L 2 ) 
is the fully discontinuous space of tensor product polynomials (shape functions in 

Qr)- 

This paper develops a second family of finite clement spaces for cubical meshes. 
This family may be viewed as an analogue of the V r A k family for cubical meshes. 
We denote the new family of elements by S r A k and define such a space for all 
dimensions n > 1, all polynomial degrees r > 1, and all form degrees < k < n. 
For 0- forms (discretization of H 1 ), S r A° is not equal to Q r , but rather to the 
serendipity elements in 2-D, their three-dimensional extension (which can be found 
in many places for small values of r and for general r in J. 2] and [8]), and to a 
recent extension to all dimensions [2]. For n-forms, S r A n uses fully discontinuous 
elements with shape functions in V r (not Q r ). 

For I-forms in 2-D, the space SrA 1 coincides with the rectangular BDM elements 
of [7]. For 2-forms in 3-D, we believe that the space S r A 2 , which has not appeared 
before as far as we know, is the correct analogue of the BDM elements on cubic 
meshes. It has the same degrees of freedom as the space given in j5] but the shape 
functions have better symmetry properties. For 1-forms in 3-D, SrA 1 is a finite 
element discretization of if (curl). To the best of our knowledge, neither the degrees 
of freedom nor the shape functions for this space have been proposed previously. 
Even for 0-forms, the spaces S r A° were only discovered very recently in higher 
dimensions. They are the generalization of the serendipity spaces given by the 
present authors in [2 . That work was motivated by the search for a finite clement 
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discretization of the de Rham complex on cubical meshes which is completed in 
this paper. 

In the remainder of the paper we will develop the S r A k spaces and their proper- 
ties in the setting of differential forms on cubical meshes of arbitrary dimensions. 
Here we explicitly describe the spaces which arise in three dimensions in traditional 
finite element terminology by giving their shape functions and degrees of freedom. 
See also Figure [T] 

The space S r A 3 . This is simply the space of piecewise polynomials of degree 
at most r with no continuity requirements. Obviously this gives a finite element 
subspace of L 2 (which, unlike the remaining spaces, is defined also for r = 0). Its 
dimension is (r + l)(r + 2)(r + 3)/6. 

The space S r A 2 . The shape functions for the space S r A 2 on the unit cube I 3 in 
three dimensions are vector polynomials of the form 

(1) U = (vi,V 2 ,V 3 ) + CUrl(lE2X3(>2 - w 3 ),x 3 xi(w 3 - Wl),XlX 2 (Wl - u> 2 )), 

with Vi,Wi G V r (I 3 ) and Wi independent of aSj. The dimension of this space is 
(r + l)(r 2 + 5r + 12)/2. As degrees of freedom for u <E S r A 2 we take 

B4 / u ■ np, p g V r (f), f a face of I 3 ; u H> / u ■ p, p e [P r _ 2 (I 3 )] 3 - 

Jf 7/3 

These degrees of freedom 

The space S r A 1 . The shape functions for the space S r A 1 in three dimensions are 
vector polynomials of the form 

(2) u = (v 1 ,v 2 ,v 3 ) + (x 2 x 3 (w 2 - w 3 ),x 3 xi(w 3 - w 1 ),x 1 x 2 (w 1 - w 2 )) + grads, 

with Vi £ V r (I 3 ), Wi G V r -i{I 3 ) independent of Xi, and s is a polynomial of I 3 
with superlinear degree at most r+ 1, where the superlinear degree of a polynomial 
is the ordinary degree ignoring variables which enter linearly (e.g., the superlinear 
degree of x 2 x 2 x 3 is 5). The dimension of this space is (r + l)(r 2 + 5r + 18)/2. As 
degrees of freedom for u £ S r A x we take 

UH J u ■ tp, p € V r (e), e an edge of I 3 ; 

J uxnp, p£ [Vr- 2 (f)} 2 , f a face of I 3 ; uh> J u-p, p6 [P r -±(I 3 )} 3 - 

The space S r A°. Finally, the space S r A° is the generalized serendipity space of 
[2]. The shape functions are all polynomials of superlinear degree at most r, and 
the degrees of freedom are the values at the vertices and the moments of degree at 
most r — 2, r — 4, and r — 6 on the edges, faces, and interior, respectively. The 
dimension of this space is 8 for r = 1, 20 for r = 2, and (r + l)(r 2 + 5r + 24)/6 for 
r > 3. 

As a consequence of the general theory below, the degrees of freedom given 
are unisolvent for all these spaces, and for any cubical decomposition T of Q, the 
assembled finite element spaces S r A k (T) have exactly the continuity required to 
belong to H l (tt), H(div, Q), i?(curl, fi), and L 2 (n), for k = 0, . . . , 3, respectively. 
In other words, S r A k (T) belongs to the domain of the exterior derivative on fe-forms. 
Moreover, exterior derivative maps S r A k (T) into S r -iA k+1 (T), so we obtain a finite 
element subcomplex 

-> S r A°(T) ^ Sr-iA^T) ^ 5 r _ 2 A 2 (T) ^> 5 r _ 3 A 3 (T) -> 
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of the de Rham complex. Finally, if we define projection operators 7r fe from smooth 
fields into the finite element spaces S r A k using the degrees of freedom, these com- 
mute with exterior differentiation. That is, the following diagram commutes: 

o — ► c°°{n) [c°°(n)} 3 [c°°(ft)] 3 c°°(n) — ► o 
*°1 ^1 " 3 1 

— * S r A°(T) ^ S r - 1 A 1 {T) 5 r _ 2 A 2 (T) ^ S,- 3 A 3 (T) — * 

where, for simplicity, we assume r > 3 (otherwise some of the spaces are undefined 
and some parts of the diagram are not applicable). 

The remainder of the paper consists of four sections. In Section 2, we recall some 
key concepts from exterior calculus, particularly the Koszul differential and Koszul 
complex of polynomial differential forms, which will be crucial to our construction. 
We introduce the concept of linear degree, and show that the subcomplex obtained 
from the Koszul complex by placing a lower bound on linear degree is exact. This 
is a key step in the unisolvence proof in Section 3. In Section 3, we define the 
spaces S r A k on an n-dimensional cube by giving the shape functions and degrees 
of freedom. We derive a number of properties from these definitions, leading to a 
formula for the dimension of the space of shape functions and a proof of unisolvence 
of the degrees of freedom. We then define projection operators w k mapping smooth 
fields into the finite element spaces 5 r A fc and show that the projections commute 
with the exterior derivative. This accomplished, the new elements fit squarely into 
the framework of the finite element exterior calculus given in [5]. Therefore the 
application of the elements to PDE problems and their numerical analysis does not 
require new ideas, and so we do not discuss that here. The proof of unisolvence in 
Section 3 hinges on the special case of functions with vanishing trace. The proof in 
that case is the topic of Section 4. 



2. Notation and preliminaries 

For n > 1, the number of dimensions, let N n — {1, . . . , n}, and let S(fc), < 
k < n, denote the set of subsets of N n consisting of k elements. For a G 
we denote by a* G S(n — k) its complement N n \ a. For p G u we write a — p 
for g \ {p} G S(fc - 1), and for q G a* we write a + q for a U {q} G S(fe + 1). If 
a C N„ and q G a* we let e(q 7 a) = (—1)' where I = #{p G a\p < q}, and set 
e(<7, p,a) := e(q, cr)e(p 1 a + q — p) for q G u* , p G u. For later reference, we note that 

(3) e(q,p,a) = -e(p,a-p)e(q,<r-p), 

which is easily verified by considering the cases q < p and q > p separately. 

We now recall some basic tools and results of exterior algebra and exterior cal- 
culus. These can be found, for example in [U Section 2]. For each a G we 
denote by a±, . . . , <7fc its elements in increasing order, and by 

dx a — dx ai A • • • A dx ak 

the corresponding basic alternator. A differential fc-form u> on a domain C K™ 
may be written as 

(4) oj = ^ u a dx a , 

o-es(fe) 
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with coefficients u) a belonging to any desired space of functions on ft, e.g., L 2 (fl). 
A 0-form is simply such a function. The exterior derivative of the differential form 
@ is 

dui = d q u> a dx q A dx a , 

trGS(fc) <?GN„ 

where d q = d/dx q . 

A differential fc-form may be contracted with a vector field on fi to give a differ- 
ential (k — l)-form (or zero if k = 0). When the vector field is simply the identity, 
the resulting operator is the Koszul differential. Equivalently, we may define the 
Koszul differential on the basic alternators by 

fe 

K,(dx ai A • ■ ■ A dx ak ) = ^Y +1 x (Ji dx ai A • ■ ■ A dx a . A ■ • • A dx ak , 

i=l 

and then extend it to a general differential form by linearity: 

fe _ 

(5) K C^2 dx„) — ui a l) i+1 ^cr i dx ai A ■ • • A dx (Ji A • • ■ A dx ak . 

(7 (7 i — 1 

The operator k is a graded differential, meaning that 

k(kuj) = 0, k(uj Arj) = kuj Arj + ( — l) k uj A K77, 

if a; is a fc-form and r\ an Worm. 

The following lemma collects formulas for kw, duj, and kcZu;. 

Lemma 2.1. 7/w is <?ii;en 6y i/ien 

(6) ™= ^ i] ( dx c , where r] C = ^2 e(q,C)x q uJ i:+q , 

Ces(fe-i) gee* 

(7) <2u; = ^ p c?a;p, w/iere v p = 2J e(e, P - q)9 q uj p - q , 

pe£(fc+i) ?ep 

(8) ftdw = 2J (J-o-dXv, where n a = ^2[x q d q uj a + '^e(q,p 1 (j)x q d p uj a+q ^ p \. 
Proof. By definition 

so 

KUJ = ^ X! 6 ^'' 7 _ q)XqU a dx a _ q . 
o-GS(fe) qeo- 

Making the change of variable £ = a — q, so 

<7GE(fc), gecr <s=>> CeS(fc-i), geC, 

we obtain (|6| . The second result is proven similarly, and the third follows from the 
first two. □ 

We now turn to differential forms with polynomial coefficients. A monomial 
in n variables is determined by a multi-index a of n nonnegative integers: x a — 
x" 1 ■ ■ ■ x" n . By a form monomial in n variables, we mean the product of a monomial 
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with a basic alternator: to — x a dx a for some multi-index a and cr c N„. The 
polynomial degree and the linear degree of m are defined 

deg m = \a\ := a,, Ideg m = i G o* \ = 1}. 

i 

Thus the linear degree of to is the degree of its polynomial coefficient counting 
only those variables which enter linearly, and excluding variables which enter the 
alternator. For ordinary monomials, i.e., 0- forms, the linear degree is equal to 
the difference between the polynomial degree and the superlinear degree which 
appeared in the introduction. 

We define W r A k — "H r A fe (IR™) to be the span of the fc-form monomials m with 
degm = r, and V r A k = V r A k {R n ) = EU-o?^ Afe to be the s P an of those with 
degm < r. If k = 0, we may simply write T-L r and V r . If £1 is a subdomain of 
E", we define V r A k (VL) to be the space of restrictions to Q, of the elements of V r A k 
(and similarly for other function spaces). Note that d maps T-L r A k into 7i r -iA k+1 
while k maps T-L r A k into H r +iA k ~ 1 . An extremely useful identity is the homotopy 
formula ([4j Theorem 3.1]): 

(9) (dK + nd)uj = (r + k)ui, uj £ U r A k (R n ). 

We extend the linear degree for form monomials to polynomial differential forms 
by defining Ideg /i for any /i € T-L r A k to be the minimum of the linear degree among 
all the monomials m in /i. We say that [i is of homogeneous linear degree equal to 
I if Ideg m = I for every monomial of fi. We denote by / K r ^A k the space of forms in 
T~L r A k of linear degree at least I. Obviously, 

(10) U nl A k (R n ) = 0, I > min(r, n - k). 

The exterior derivative d may decrease the linear degree of a polynomial differ- 
ential form, but k and ndn do not. 

Lemma 2.2. For any lj £ V r A k , ldegfaj > ldegw and Idegndnuj > ldegw. 

Proof. If m is monomial of u> and / = Ideg to, then it follows directly from the 
definition ([5| that the monomials of nm are of linear degree I and/or I + 1, so 
IdegKTO > I > Ideg a;. Since every monomial of kuj is a monomial of Km for some 
monomial m of lj, this implies the first inequality. For the second we use the 
differential property of n and the homotopy formula to see that ndnm = {nd + 
dn)nm is a multiple of nm. Therefore Ideg ndnm = \degnm > ldegcj, which gives 
the second inequality. □ 



In view of Lemma |2.2| for each I > 0, we obtain a complex 

(11) • • • A H r - U A k+1 A H r jA k A ••• . 

When I — this is the Koszul complex, and exactness follows from the homotopy 



formula. In fact, the complex (11 1 is exact for all I > 0. 
Theorem 2.3. For r >1, < I < r, < k < n, the sequence 

n r ^ u A k+1 A n rJ A k A n r+ uA k ^ 

is exact. Equivalently, dim K( - H r -i ! ;A fc+1 ) + dim K(T-l r jA k ) = dim'H r) ;A fc . 

The proof is due to Scot Adams and Victor Reiner [T]. Its main ingredient is 
contained in the following lemma. 
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Lemma 2.4. Let r > \, I > 1, and < k < n. Suppose that p £ 'H r A k is of 
linear degree at least I — 1 and Kfi is of linear degree at least I. Then there exists 
v £ H r _iA fc+1 such that p — kv is of linear degree at least I. Further, if p £ % r A n 
is nonzero, then Kp is of linear degree 0. 

Proof. For the final statement, concerning n-forms, we write p = pdx\ A • • • A dx n 
where p £ H r - Since r > 1, p is not constant. But then Kp = l) i+1 pxi dx\ A 

• • • A dxi A • • • A dx n is easily seen to be of linear degree 0. 

For < k < n, the proof hinges on a canonical form for an element of p £ T-L r K k 
which we establish before proceeding. Let us say that a form monomial x a dx a is 
full if a C supp(a), the support of a. To each of the monomials m — x a dx a of 
p £ TirA k , we associate the increasing sequences p and r with p — crnsupp(a) and 
t = a\ supp(a). Then m = ±rj A dx T where r\ = x a dx p is a full form monomial 
which is independent of the r variables (that is, supp(a) U p is disjoint from r). 
Note that ldegry = ldegrn. Finally, in the expansion of p as a linear combination 
of its monomials, we gather together the terms with the same r = a \ supp(a), and 
in this way write 

(12) p = Tjr A dx T , 

tCN„ 

where r\ T £ T-L r k k ~# T is independent of the r variables, and has all of its monomials 



full. The expression on the right-hand side of (12 1 is the desired canonical form of 
p. 

Now we proceed with the proof of the lemma. We consider first the special case 
in which p is of homogeneous linear degree I — 1, and, in this special case, we use 
induction on I, the case / = being known (exactness of the Koszul complex). 
Expressing p in the canonical form ( [l2| , we have ldeg77 T = I — 1 for each r (for 
which the coefficient r\ r does not vanish). Now 

(13) up = kt/ t A dx T + ±rj T A ndx T . 

The first sum is of homogeneous linear degree I — 1 and the second of homogeneous 
linear degree I. Since we assumed that up is of linear degree at least I, the first 
sum must vanish. But this sum is in canonical form, so we conclude that nr) T = 
for each r. Invoking the inductive hypothesis, we can write r\ T — nv T where v T has 
linear degree at least I — 1. Let 



" = £■ 

Then 



Adx T £ H r -iA k+1 . 



p — KV 



r] T A dx T — kv t A dx T — ±v T A Kdx T — — ±v T A Kdx 7 



Expanding the right-hand side into monomials, we see that each has ldeg at least I, 
so this completes the proof under the assumption that p is of homogeneous linear 
degree I — 1. 

Next we turn to the general case, in which p is of linear degree at least / — 1. We 
may split p as p! + p" with p! of homogeneous linear degree I — 1 and p" of linear 
degree at least I. Then Kp! splits into a part of homogeneous linear degree I — 1 and 
a part of homogeneous linear degree /, while ldeg(K/i") > I. Since, by assumption, 
\deg(Kp) > I, the part of Kp' with linear degree equal to I — 1 must vanish. That is, 
Ideg(Kp') > I. Therefore, we may apply the result of preceding special case to p! to 
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obtain v € 'H r -iA fe+1 such that ldeg(// — kv) > Z. Then fi — kv = (fjf — kv) + fi" 
is of linear degree at least I. □ 

Finally, we give the proof of Theorem |2.3| 

Proof. The result is certainly true for I = 0, so we may assume 1 < I < r (and so 
r > 2). Suppose co g H r ,iA k with kw = 0. We must show that there exists r\ € 
% r -iA fc+1 with linear degree at least I, such that kt/ = cj. Now /i := dui/(r + k) 6 
■H r -iA fc+1 is of linear degree at least i — 1 and satisfies K[i — (Kd + dn)ui/(r + k) = uj 
by (J9|. If k = n — 1, the final sentence of the lemma insures that /i = 0. For 
0<fc<rt~ 1, we apply the lemma with r and fc replaced by r — 1 and fc + 1, 
respectively, and conclude that there exists v € H r -2^- k+2 such that i] := /i — kis is 
of linear degree at least Z. Clearly k?7 = k/i = cj. □ 

3. The S r A k spaces 

In this, the main section of the paper, we define the polynomial spaces S r A k (M. n ) 
we shall use as shape functions (see p"7| ) below) and the degrees of freedom for 
these (see (pTJ). We derive a number of properties of these polynomial spaces in 
Theorems |3.2| through |3.5| and use them to verify unisolvence in Theorem |3.6| 

The space of shape functions will consists of polynomials of a given degree plus 
certain additional terms of higher degree which will be defined in terms of the 
following auxilliary space: 

J r A k (R n ) = J2^n r+ i- 1 . l A k+1 (R n ). 
i>i 



In view of ( 10 ), the sum is finite and 



(14) J r A k (R n ) c n.+„_ fc _ 1 A fc (r i ). 

Moreover, the sum is direct, since the polynomial degrees of the summands differ. 
The following proposition, which follows directly from the definitions, helps to 
clarify the meaning of this space. 

Proposition 3.1. 1. The space Yli>i Hr+i-i,iA k (M. n ) is the span of all k-form 
monomials m with deg m > r and deg m — ldeg m < r — 1 . 

2. The space J r A (R") is the span of Km for all (k + l)-form monomials m with 
deg m > r and deg m — ldeg m < r — 1 . 



For several values of k, this space can be described more explicitly. By (10 1, 
(15) J r A k {R n ) =0, for k = n or n - 1, 

while 

J r A n - 2 (R n ) = k-H^A"- 1 ^"). 
Now, 7i r ,iA n ^ 1 (R n ) is the span of the monomials XiWi6i, where Wi £ T-L r -i(R n ) is 
independent of Xi and ftj := (— \) l ~ 1 dxi A ... A dx % A ... A dx n . We then have 

j<i j>i 

with 

9 i:j = {-l) i+j dx 1 A...Adx { A...Adxi A... A dx n . 
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Therefore, 

(16) J r A n - 2 (R n ) 



{ XiXj(wi — Wj)8i,j | Wi E "H r -i(]R n ) independent of Xi }. 



l< 3 



Finally we identify J r A°(R n ). By Theorem [2^3 

= J2l>l Hr+l,lA°(M. 



in the case k = 0, we see that 
this space is the span of 
ldeg, is at most 



3.1 



J r A u (R n ) = Ei>i^r+i,iA u (IR"). By Proposition 
monomials of degree > r whose superlinear degree, that is, deg 
r. 

We can now define the space of polynomial fc-forms which we use for shape 
functions: 



(17) S r A (R n ) := V r A*(R n ) + J r A (R n ) + dJ r+1 A k 

defined for all r > 1 and all < k < n. From JTil 



k-it 



(18) 



S r A k (I 



C V r 



Note that, in case k = 0, the final term in (17) vanishes, and, by the charac- 
terization of J r r A°(M n ) just derived, S r A°(R n ) consists precisely of the span of all 
monomials of superlinear degree at most r. This is exactly the serendipity space as 
defined in [2] . In this case, ( 14 ) gives the sharper degree bound 

c P r +«-iA°(I 



(19) S r A°(l 

Another case in which the expression for S r A k (I 
k = n. By (pi, 



l n ) can be simplified is when 



S r A n (R n ) 
When k = n — 1, we have, by ( p| , 

5 r A n_1 (R") :=V r A n - 



V r A n (l 



L (R n ) + dJ- r+ iA n - 2 (R n ), 



where the last space is characterized in In the case of three dimensions, this is 
the formula given in the introduction, stated in the language of exterior calculus. 
In a similar way, we recover the formula ^ for the 3-D if (curl) elements discussed 
in the introduction. 

We now derive several properties of these polynomial spaces. The first limits the 
monomials that appear in the polynomials in 5 r A fe (]R"). 

Theorem 3.2 (Degree property). For any n, r > 1 and < k < n, the space 
S r A k (JH ri ) is contained in the span of the k-form monomials m of degree at most 
r + n — k — Sko for which 

(20) degm — ldeg to < r + 1 — <5fco- 

Proof. The bound r + n — k — Sko on the degree is given in ( 18 ) for k > and 
in (19 1 for k = 0, so we need only show (20). If to is a monomial of an element 
of V r A (R"), then degm < r and ldegm > 0, so (20) holds with r on the right- 
hand side. If to is a monomial of an element of t 7 r A' c (R"), then to occurs in the 
expansion of up, where p is a (k + l)-form monomial with degp — ldegp < r — 1 
(Proposition 3.1). Then degm = degp+ 1 and, by Lemma 2.2 IdegTO > ldegp, so 
again degm — IdegTO < r. Finally, if k > and to is a monomial of an element of 
e? l 7 r +iA fe ~ 1 (K n ), then by the argument just given, to is a monomial of dq where q 
is a (k — l)-form monomial with degq — ldegq < r + 1. Since degm = dcgq — 1 
and ldegm > ldcgq — 1, we get (20). □ 
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A crucial property of these polynomial form spaces, is that they can be combined 
to form a subcomplex of the de Rham complex. 

Theorem 3.3 (Subcomplex property). Let n,r > 1, and let < k < n. Then 

dS r+ iA fe - 1 (R") c S r A k {R n ). 



Proof. With reference to ([Tr}, we note that dV r+ iA k - 1 (R n ) C V r A k (R n ), and 
d(dJ r+2 A k (R n )) vanishes, so it suffices to prove that dJ r+ iA fc - 1 (M") C S r A k (R n ) 1 
which is immediate from (17). □ 

We next observe that the spaces increase with increasing polynomial degree. 

Theorem 3.4 (Inclusion property). Let n, r > 1, and let < k < n. Then 

S r A k (R n ) cS r+1 A k (R n ). 



Proof. We must show each of the three summands on the right-hand side of (17) 
is included in S r+1 A k (R n ). Clearly 

V r A k (R n ) c V r+1 A k (R n ) c S r+1 A k (R n ), 

which establishes the first inclusion. Next, we show that 

J r A k (R n ) c V r +iA k (R n )+J r+1 A k {R n ) c S r+1 A k (R n ). 

By Proposition 3.1 elements of J' r A k (R n ) are of form Km with m a (k + l)-form 



monomial with degm > r, degm — ldeg to < r — 1. By the homotopy formula 
nm is a constant multiple of up with p = dnm. We have degp = deg to and ldegp > 
ldeg to — 1. Then degp — ldegp < deg m — ldeg m + 1 < r. If dcgp = r, deg Km — 
r + 1 an d Km e T' r +iA k (R n ). On the other hand, if degp > r + 1, by Proposi- 



tion 3.1 Km G J~ r +iA (R n ). This establishes the second inclusion. To complete 
the proof, we show that dJ r+1 A k - 1 (R n ) C S r+1 A k (R n ). Since J r+1 A k ^ 1 (R n ) C 
iS r +2A fe_1 (IR™) (by the inclusion just established), we infer from the subcomplex 
property that dJ r+1 A k - 1 (R n ) C dS r+2 A k - 1 (R n ) C 5 r+ iA fe (IR™). □ 

The third property of the S r A k spaces we establish concerns traces on hyper- 
planes. Consider a hyperplane / of R n of the form x.- L = c for some 1 < i < n 
and some constant c. The variables Xj, j ^ i, form a coordinate system for /, 
so we may identify / with and consider the space S r A k (f). It is a space of 

polynomial fc-forms on /, and so vanishes if k = n. Next, we consider the trace on 
/ of a differential form in n variables (defined as the pullback of the form through 
the inclusion map / R n ). Let a £ £(£;), and let oj a be a function of n variables. 
Then 

0, i € cr, 

(tr/ cjcr) dx a , i ^ (7. 

In the last expression, tr j u a denotes the function of n — 1 variables obtained by 
setting Xi — c and we view dx a as a basic alternator in the n—1 variables Xj, j =/= i. 
The trace property states that if u £ S r A k (R n ), then tifU, which is a polynomial 
fc-form on /, belongs to S r A k (f). 

Theorem 3.5 (Trace property). Let n, r > 1, < k < n, and let f be a hyperplane 
ofR n obtained by fixing one coordinate. Then 

tr f S r A k (R n ) C S r A k (f). 

(This inclusion will be shown to be an equality in (27) below.) 
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Proof. Without loss of generality, we assume that / = { x € R n \x% = c}. First 
let us comment on the Koszul operator applied to a polynomial differential form 
on /. Such a form may be written as a linear combination of monomials x a dx a 
where a\ = and 1 ^ a. Referring to ^ we see that, if we view x a dx a as a form 
monomial in n variables and take the Koszul differential, the result is the same as 
if we view it as form monomial in n — 1 variables and take the Koszul differential. 
Thus we need not distinguish between the Koszul differential on K™ and that on /. 

We will prove the theorem by induction on k. For k = we recall that S r A°(W l ) 
is the serendipity space spanned by the monomials of superlinear degree at most 
r, and, of course, the superlinear degree does not increase when taking the trace. 
Hence, tr/5 r A°(K") C 5 r A°(/). 

To prove the theorem for k > 0, assume that it holds with k replaced by k — 1. In 
light of (JTYJ) we need to show that the traces of each of the three spaces V r A k (H n ), 
J r A k {R n ), and d^+iA*- 1 ^") are contained in 

S r A k (f) = V r A k (f) + J r A k (f) + dJ r+1 A k -\f). 

For the V r A k (R n ), this is evident, since try-p r A fc (R n ) C V r A k (f). 
Next, we establish that try dJ r+1 A k ~ 1 (R n ) C S r A k {f). Indeed, 

try dJ r+ iA k - 1 (R n ) = dtx } J r+l A k - l ^ n ) 

C dtr f S r+1 A k -\R n ) c dS r+1 A k -\f) c S r A k (f), 
where we have used, in turn, the commutativity of the trace with exterior differen- 



tiation, (17), the inductive hypothesis, and Theorem 3.3 



By Proposition 3.1 in order to show that try J' r A k (W l ) C S r A k (f), and so 
complete the proof, it suffices to show that try Km E S r A k (f) whenever m is a 
(k + I)-form monomial with degm — ldegm < r — 1. We write m as x a dx a , and 
consider separately the cases 1 ^ a and 1 6 a. 

Assuming 1 ^ <j, let p be the (fc + l)-form monomial obtained restricting m to 
/, i.e., by setting x± — c in the coefficient x a . Then try Km — up. If m is linear 
in xi, then degp = degm — 1 and ldegp = ldegm — 1. Otherwise degp < degm 
and Idegp = ldegm. In either event, degp — ldegp < degm — ldegm < r — 1, so 
up G S r A k (f) (again using Proposition 3.1 ). 



Assuming, instead, that 1 € cr, we may write 

m = x a dx a = x^x^ dx\ A dx T , 

where /? is a multi-index with f3\ = and r € E(fc) has t± > 1. Then 

p := try nm — c° 1+1 x^ dx T 

is a /c-form monomial independent of x\ with degp < degm and ldegp = ldegm. 
We are trying to show that p e S r A k (f). This is obvious if degp < r, so we may 
assume that degp > r. We shall show that both Kdp and dKp belong to S r A k (f), 
which suffices by ([9]). 

Now deg dp = deg p — I < deg m — 1 and ldeg dp > ldeg p — I = ldcg m — 1 , so 
deg dp — ldeg dp < degm — ldegm < r — 1. Therefore Kdp £ J r A k (f) C S r A k (f), 
as required. 

Finally, we show that Kp e j7 r +iA fe ~ 1 (/), whence d«p G S r A k (f) as well. By 
Proposition |3.I| this holds, since degp > r+l and degp— ldegp < degm— ldegm < r 
(even < r — 1). This concludes the proof. □ 
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Having denned the space S r A k (M. n ) of polynomial differential forms, we turn now 
to the definition of the associated finite element space on a cubical mesh. As usual 
the finite element space is defined element by element, by specifying a space of 
shape functions and a set of degrees of freedom on each cube T in the mesh. (More 
generally the element T may be a right rectangular prism, that is the Cartesian 
product of n closed intervals of positive finite length.) As shape functions on T 
we use S r A k (T), the restriction of the above polynomial space to the cube. The 
degrees of freedom have a very simple expression. Writing A<j(T) for the set of 
d-dimensional faces of T, they are given by 
(21) 

H H- f tifn A v, v G V r - 2(d -k)A d ~ k (f), f G A d {T), k < d < min(n, [r/2\ + k). 
Jf 

Since 

di^^A-'tf) = dlmP^CA * („ t ,) - ('• - 7 2k ) ( d k ) 

for / a face of dimension d, and since there are 2™~ d Q) d-dimensional faces of an 
n-cube, the number of degrees of freedom in (21) is given by 

min(n,Lr/2j+fe) / \ / . , o; \ /, 

We now turn to one of the main results of this paper, the proof that the degrees 
of freedom (21) are unisolvent for S r K k (T). 

Theorem 3.6 (Unisolvence) . Let n, r > 1 and < k < n, and let T be a cube in 
W l . Then 

(1) dim5 r A fc (T) is given by ((22). 

(2) If /i G S r A k (T) and all the degrees of freedom in (21) vanish, then /i = 0. 

Using the trace property, we will reduce the proof of (2) to the case where /j, 
belongs to the space 

S r A k (T) :={/ie S r A k {T) \ tr/ fi = on each face / G A n _ x (T) }, 

the subspace with vanishing traces. This case is given in the following proposition. 

Proposition 3.7. If fi G S r A k (T) and 

(23) / nAiy = 0, ^e7> r - 2 („- fc) A"- fc (T), 

Jt 

then /i vanishes. 

We defer the proof of this proposition to Section |4j Now, assuming this result, 
we prove Theorem |3.6| 



Proof. We begin with the first statement of the theorem. Since J r A k (T) is in the 
range of k, and since the homotopy formula implies that no nonzero differential 
form is in the range of both k and d, the sum on the left of (17) is direct. The 
homotopy formula implies as well that d is injective on the range of n. Therefore 

dim5. r A fe (r) = dimV r A k (T) + ^dimKU r+ i- 1 jA k+1 (T) + J2dimKH r+ i t iA k (T). 

/>! 1>1 
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Applying Theorem |2.3[ this becomes 



dimS r A fc (T) = dim-p,.A fe (T) + ^ dimH r+M A fc (T) 

= dim[P r A fc (T) + ]TH r+M A fe (T)]. 



i>i 



i>i 



The space in brackets is exactly the span of the fc-form monomials m with deg m - 
ldegm < r, and hence we need only count these monomials. This gives 



dimS r A k (T)= r)#A(r,k,n) 



(24) 



where is the number of basic alternators and A(r, n, k) is the set of monomials 
p in n variables which are linear in some number I of the first n — k variables, 
with degp — I < r. Now we count the elements of A(r, n, k). For any 
monomial p in n variables let J C N„_fe be the set of indices for which Xi enters p 
superlinearly, let d > be the cardinality of J, and for i € N„_fc \ J, let a, = or 
1 according to whether p is of degree or 1 in Xj . Then 

a**** n 

where g is a monomial in the d variables indexed by J and the last k variables. With 
I the number of aj equal to 1, we have degp = 2c? + deg q + l. Thus degp — I < r if 
and only if deg q < r — 2d. Thus we may uniquely specify an element of A(r, k, n) 
by choosing d > 0, choosing the set J consisting of d of the n — fc variables (for 
which there are ( n ^ fe ) possibilities), choosing the monomial q of degree at most 
r — 2d in the d + k variables ( ( r ^ 7.^ ) possibilities) , and choosing the exponent aj 



P = 



(n 



to be either or 1 for the n — fc 



-k-d 



(25) 



d remaining indices (2 

min(n— fc, L r /2j ) 

#A(r,k,n)= 2 ™ 

ci=0 

min(n,Lr/2J+fc) 

= E 2 " 



i, — k — 



possibilities). Thus 



n — k 
d 



n — k 
d-k 



r — d + k 
d + k 

-d + 2k 



where the second sum comes from a change of the summation index (d — > d — k). 



Substituting ( 25 ) into ( 24 1 and using the binomial identity 

'n\ fn-k\_ fd\ fn 
Ml \d-kl ~ [kl \d 



we conclude that 



(26) 



dim<S r A fc (T) 



min(n,Lr/2j+fe) 

E 2 "" 



d=k 



d + 2k 



This completes the proof of the dimension formula for 6> r A fc (]R"). 

The proof of unisolvence is easily completed based on the trace property and 
Proposition |3.7| We use induction on the dimension n, the one-dimensional case 
being trivial. Suppose [i £ S r A k (T) and all its degrees of freedom vanish. For any 
face / of dimension n — 1, try /i £ S r A k (f) and all the degrees of freedom for it 
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vanish. By induction try /i = on /. This implies that /i € S r A k (T), and we invoke 
Proposition |3.7| to conclude that /i vanishes identically. □ 



We remark that, as a corollary of unisolvence, we may strengthen the result of 



Theorem 3.5 to equality: 

(27) tr/5 r A fc (M n ) =S r A k (f). 

Indeed, let T be a cube with one face contained in the hyperplane /. Then T D f 
is an {n — l)-dimensional cube and any v £ S r A k (f) is uniquely determined by 
the degrees of freedom for the S r A k (T n /). Now we may determine an element 
jU G S r A k (R n ) by assigning the degrees of freedom for the space S r A k (T) arbitrarily. 
In particular, we may choose the values of those degrees of freedom associated to 
the face Tf]f and its subfaces to be the same as those for v. Then tr/ fi € S r A k (f) 



(by Theorem 3.5), and try /i and v have identical degrees of freedom, and so they 



are equal (by Theorem 3.6) 



With the definition of the spaces complete, we use the subcomplex property, 
Theorem |3.3[ to dehne a subcomplex of the de Rham complex on the cube: 

R S r A°(T) -A Sr-iA^T) ... -i* 5 r _„A"(T) — ► 0. 

To show that this complex is exact, we define the canonical projection 

n k : CA k (T) ->5 r A fc (T), 

associated to the unisolvent degrees of freedom. That is, Tr k /i g S r A k (T) is deter- 
mined by the equations 

i tXf{-K k ix)hv = Jtr fl jLhv, ^P^^fc)^^/), f€A d (T), 

k <d< min(n, \r/2\ + k). 

Then, the following diagram commutes: 

R C°°A (r) C°°A 1 (T) ••• C°°A n {T) — > 



R 5 r A°(T) -A Sr-iA^T) -1+ ... -1> 5 r _„A"(T) — > 

The proof of commutativity is based on two basic properties of differential forms: 

(1) the commutativity of trace and exterior differentiation, try do; = dftifuj, and 

(2) integration by parts, which for differential forms can be written 

duj A r? = (-l) fe_1 / LuAdr]+ / tr 9n uj AtrgnT], 
in J a, Jan 

for a fc-form u) and an (n — k — l)-form rj on an rt-dimensional domain fi. See [U 
Lemma 4.24] for the same argument applied to simplicial elements. Since the top 
row of the diagram, the de Rham complex on the cube, is exact, the commutativity 
of the diagram implies that the bottom row is exact as well. 

Having defined the finite element space S r A k (T) on a single cube T and estab- 
lished its properties, the space S r A k (Th) associated to a cubical mesh Tk is defined 
through the usual finite element assembly. In view of the unisolvence result Theo- 
rem 



3.6 and the trace result ([27]), the degrees of freedom associated to a face of the 
cube and its subfaces determine the trace of the finite element differential form on 
the face. It follows that S r A k {T h ) C HA k (fl) (see H Section 5.1]). 
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4. UNISOLVENCE OVER THE SPACE WITH VANISHING TRACES 



We conclude the paper with the proof of Proposition |3.7| which is based on the 
following lemma. 

Lemma 4.1. Suppose that rj = 53<reS(fe) d%a where, for each a G X]„ is 

a homogeneous polynomial which is superlinear in all the a* variables. Further 
suppose that ldeg nn > 1 and ldeg ndrj > 1. Then r\ = 0. 

Proof of Lemma \4J\ By ([8| we have ndn = Ylo-eEtk) dx&, wnere 

(28) \l a = ^ [XqdqVo + ^2 P ' CJ ) X q d P 1 la+q-p\ ■ 

Since ldeg ndn > 1, each monomial of the polynomial \i a is linear in at least one a* 
variable. 

For any a G S(fc) and any subset r of a, let 5*^ be the span of the (ordinary, 
0-form) monomials which are independent of the r variables but depend on all the 
other a variables, and let S* denote the span of the monomials which are superlinear 
in all the a* variables. Denote by PJ the projection onto S£. That is, if p = J2Pm m 
where the sum is over all monomials m and the coefficients p m are real numbers, 
all but finitely many zero, then P£p := X^meS T Pm TO - Similarly denote by Q T a the 
projection onto S^Pi S*. We now calculate the result of applying Q T a to both sides 



of ( 28 ) . First we note that 



(29) Q T a ^ = 0, 

since every monomial of /i a is linear in at least one a* variable, and so none of them 
belong to S*. Next, for each q G a*, 

(30) QKxqdqTiv) = PZ{x q d q T) a ) = x q dq(PZr)o), 

with the first inequality holding since each monomial of r) a , and so of x q d q r\ a is 
superlinear in all the a variables. Now we determine the action of Q T a on the terms 



of the second sum on the right-hand side of (28 1. For any q G a* and p G er, we 
claim that 



(31) Q T a (x q d p r) a+q - p ) 



0, per, 
^x q d p (PZ +q _ p r)a+q-p), p€o-\r. 

Indeed, r] a+q _ p is superlinear in x p , so every monomial of x q d p r\ CT j rq - p depends on 
x p . This implies that the projection is in the case p G r. In case p G a \ r, we 
write x q d p ri a+q - p = ^ x q d p m, where the sum is over the monomials m of r\ a + q - p . 
Since neither p or q belongs to r, the monomial m' := x q d p m is independent of 
the r variables if and only if same is true of the monomial m, and, since m! always 
depends on x p , it depends on all the a \ r variables if and only if m depends on all 
the (a — p)\r variables. Further, m! is always superlinear in all the a* variables 
except possibly x q , and it is superlinear in x q if and only if m depends on x q . In 
short, m' belongs to S£ n S* if and only if m belongs to S^ +q _ p . This completes 



(301, and Oik, that 



the verification of (31). Thus the application of Q T a to (28) gives, in light of (29), 



(32) ~J2[x q d q (P^r] a )+ e (<l,P,v) X qdp( P *+q- P V<J+q~p)] = 0. 

qCv* pGcr\r 
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We now claim that, for any a £ and r C u, that 

—pVcr+q—p) — ( c + 
q£<y' pEcr\T ieer\r 

where c = #(cr \ r). Assuming this, we have from ( [32] ), 

= (c+ Y + E X A)(KV*) = (c+ E ^)(P;r? CT ) = (c + r + l)(F ff V), 

i£a\r q£tr* iET* 

where we have used Euler's formula for homogeneous polynomials (i.e., the ho- 
motopy formula for 0-forms) in the last step. Thus P^,r\ a vanishes, and so r}„ = 
^ tCo . P^riu also vanishes. Since a £ S(fe) was arbitrary, this implies that r\ = 0. 
Thus it remains only to prove ( 33 ) . By the first formula of Lemma |4j n-q — 

Eces(fc-i) w C dx C where 

W C = E e (l-'OXqVC+q- 

Now let £ 6 E(fc — 1) and t C (, and denote by i?£ the projection onto SJ- n S£ 
(the span of monomials which depend on all the £ variables except the r variables 
and which arc super linear in the £* variables). By the hypothesis that ldegK?7 > 1, 

RIloq = 0. 

Next we compute R^(x q rj^ +q ) for q £ (* . If m is a monomial of rj^+ q , then m' = 
x q riQ +q depends on all the £ variables except for the r variables if and only if the 
same is true of m. Moreover, m! is superlinear in all the Q* variables if and only 
if m depends on x q (since m is superlinear in all the £* variables with the possible 
exception of x q ). Thus 

Rl( x qVC+q) = X q P C+qVC+1 

Combining the last three displayed equations, we obtain 

Now choose some p £ C* an d differentiate this equation with respect to x p to get 

efe 9 p( X P P (+P^+p) + E £ ( ? ' X qdp( P C+qVC+q) = °> 

or, after rearranging, 

(34) P( +p VC+p + x pdp(Pc+ p V<:+p) = - E e (^C)e(<7,C)^?5 p (P c r +9 ?7 C+g ). 

<?eC*-p 

For any cr £ S(fc) and any p £ a set £ = er — pe E(fc — 1). Then we can rewrite 
(34 1 in terms of a: 

KVa + x pdp(PvVa) = - E ° ~ P ^ 6 ( q > ° ~ P) X 1 d p( P v+q-p r l°+q~p)- 

Finally, taking any t £ a £ we sum over p G cr \ t to obtain 

C 0<x + E X P d p(KVa) = - E E e (P^-P) e (9^-p) X q9p(P^ +q -pVcT+q-p)- 

where c = #(cr\r). In light of Q, this establishes (33), and so completes the proof 
of the lemma. □ 
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Finally we give the proof of Proposition |3.7[ 

Proof. By dilating and translating, it suffices to prove the result when T = I n with 
I = [-1,1]. Let 

crGS(fe) 

be a polynomial differential form on the cube I n . Then tr f fi vanishes on the faces 
X{ = ±1 if and only if for each a such that i (£ a, 1 — xf divides \i a . Now suppose 
that n € S r A k (I n ), so that tr/ /i vanishes on all the faces of the cube. It follows 
that 

iiu = fj-a ]J ( X ~ X i ) 

for some polynomial p, a . The monomial expansion of fj, then contains the form 
monomial ™ CT IIjeer* x i dx a , where m a is any monomial of highest degree of fi a 



The linear degree of this form monomial is 0, so, by the degree property (201, its 
degree is at most r + 1. 

Having established that fj, is of degree at most r+1, let r\ be its homogeneous part 
of degree r + 1. We have ldeg?7 = 0. Now we may match terms in the definition 



( |T7| of S r A k (n) to obtain that i] = kv + oIklo, for some v G / H r ^A k+1 (I n ) and 



£ H r +i,iA k {I n ). Therefore, ldegK?7 = ldcg ndnuj > ldeg u) > 1 and ldeg ndr/ = 
\degndnv > ldegf > 1, where we have used Lemma |2.2| 

By Lemma |4.1| r\ = 0, hence the monomial of highest order in the expansion 
of /i CT is of degree at most r. It follows that jl a is of degree r — 2(n — k) for each 
a e S(fc). We can then choose the test function v = J2aeY;(k)(~ l) sgn ( <T ' <7 ^jl a dx a * 



in (23) to conclude that /i vanishes. □ 
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